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Abstract 

Pixel lensing, gravitational microlensing of unresolved stars, is potentially much 
more sensitive and much more widely applicable than is generally recognized. I 
give explicit expressions for the pixel noise induced by a time- variable PSF, by 
photometric and geometric misalignment, and by discrete pixelization, and I show 
that these can all be reduced below the photon noise. Pixel lensing can be divided 
into two regimes. In the "semi-classical" regime, it is similar to classical lensing 
in that it measures the time scale of individual events. In the "spike" regime, it 
measures the total optical depth but not individual time scales. I present simple 
expressions for the boundary between the two regimes and for the event rate in 
the latter one. These expressions can be used to quickly classify all potential 
pixel lensing experiments. Pixel lensing can measure the luminosity function as 
well as the mass function of stars in target galaxies to a distance of a few Mpc. 
Future space-based pixel lensing could be ~ 5 times more effective in the infrared 
than in the optical, depending on developments in detector technology. Pixel 
techniques can also be applied to non-pixel-lensing problems including the search 
for unresolved variable stars and follow up observations of lensing events found in 
classical lensing searches. To benefit fully from pixel-lensing techniques, follow-up 
observations should have resolutions of at least 5 pixels per FWHM. 

Subject Headings: gravitational lensing - stars: masses - stars: luminosities - 
techniques: photometric 



1 Alfred P. Sloan Foundation Fellow 



1 



1. Introduction 



Pixel lensing is gravitational microlensing of unresolved stars. In classical lens- 
ing surveys of the type proposed by Paczyhski (1986), one monitors large numbers 
of resolved stars in hopes of finding a few which are microlensed by foreground 
objects. The magnification A[x) depends only on the ratio x = 6 /Be of the angular 
separation of the source and the lens to the Einstein radius, 
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M is the mass of the lens, and JqI, Jig, and Jqs are the distances between the 
observer, source, and lens. At least four such surveys are currently being carried 
out toward the Large and Small Magellanic Clouds (LMC and SMC) and the 
Galactic bulge and at least 100 candidate microlensing events have been reported 
to date (Alcock et al. 1993, 1995a, b; Aubourg et al. 1993; Udalski et al. 1994a; 
AUard 1995). The motion of the lens relative to the observer-source line of sight 
is usually well-approximated as a uniform velocity, v, in which case the separation 
is a function of only three parameters, 

where to is the time of maximum magnification, fi is the impact parameter in units 
of ^e, and LO~^ is the Einstein crossing time. The optical depth r can then be found 
directly from the observed time scales '-o~^ ^ 
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where N is the total number of stars observed, T is the duration of the observations, 
and ej are the efficiencies of detection. 

Note that the event rate Fo per observed star is a function not only of r but 
also of the typical time scale, 

Fo = -LOT. (1.5) 

TT 

While classical lensing searches are a fantastic success, they are fundamentally 
limited to galaxies with substantial numbers of resolved stars, namely the three 
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galaxies currently being searched. There is a great deal that one could hope to 
learn from lensing surveys of other galaxies. For example, the halo of M31 could be 
probed for lenses by monitoring stars in the disk of M31 (Crotts 1992; Baillon et 
al. 1993). Because of the M31's high inclination, the lensing rate should be much 
higher toward the far side of the disk compared to the near side (Crotts 1992), 
which would be an unambiguous signature of lensing. One could measure the mass 
function of the bulge of M31 by monitoring the stars in the M31 bulge. One could 
look for intra-cluster lenses in the Virgo cluster by monitoring M87 (Gould 1995b), 
or even study the star-formation history of the universe by searching for lenses at 
cosmological distances (Gould 1995a). All of these projects require pixel lensing, 
the detection of lensing of unresolved stars. 

In this paper, I present a general theory of pixel lensing. I show that there 
are two distinct regimes of pixel lensing which I label "semi-classical" and "spike". 
The semi-classical regime is similar to classical lensing. Indeed for nearby galaxies 
like the LMC, semi-classical pixel lensing may be regarded as a simple extension 
of classical lensing beyond the so-called "crowding limit". Moderately distant 
galaxies like M31 are generally still in the semi-classical regime. Classical lensing 
in these galaxies is virtually impossible, but pixel lensing yields very much the 
same type of information as classical lensing does in galaxies that are nearer. In 
particular it is possible to measure the individual time scales of the events and so 
use equation (1.4) to determine the optical depth. It is also possible to measure the 
unlensed fluxes of the source stars and thus the luminosity function (LF). In more 
distant galaxies, pixel lensing enters the spike regime. It is no longer possible to 
measure the time scales of individual events nor the unlensed fluxes of the sources. 
Remarkably, however, it is still possible to measure the optical depth. 

I present simple formulae for the boundary between the two regimes and for the 
number of events which can be detected in the spike regime with a given signal to 
noise, assuming that the noise is dominated by photon statistics. There are several 
sources of systematic error which could in principle take precedence over the photon 
noise. These include noise introduced by variation in the seeing, by problems in 
aligning the images geometrically and photometrically, and noise from discrete 
pixelization of the point spread function (PSF). I derive explicit formulae which 
allow one to determine whether these systematic sources of noise are significant 
relative to photon noise. 

The formalism presented here should prove useful in optimizing the design 
of pixel lensing experiments, in comparing different possible experiments, and in 
developing new pixel lensing ideas. I give a few illustrative examples of these 
applications. 

Two groups are currently carrying out pixel lensing searches toward M31, 
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Crotts and Tomaney (Crotts 1992) and AGAPE (Baillon et al. 1993, Melchior 
1995). Both groups have for their initial search adopted a "threshold" approach 
to identifying events. In this approach, one demands that the event rises above 
the background by a minimum number of a for several consecutive measurements. 
CoUey (1995) has analyzed the detectability of events toward M31 by simulating 
this threshold approach. 

The framework adopted here is quite different. I analyze detectability from 
the standpoint of the total signal to noise of the event assuming that photon 
statistics are the only source of noise and that the entire sequence of images can be 
searched with a set of optimal event filters. Many will undoubtedly regard these 
assumptions as symptoms of the wild optimism and hopeless naivete of a theorist. 
I was led to this approach after visiting the AGAPE group and being shown the 
difference between two deep images of M31 taken in similar seeing. The individual 
images were highly mottled with surface brightness fiuctuations, but the difference 
contained only a single "star" plus photon noise. Subsequently, Tomaney & Crotts 
(1996) have used convolution to produce differences of images taken in substantially 
different seeing. These difference images are also near the photon-noise limit. In 
one, for example, the empirically calibrated fiux error is only R = 24.6 in a region 
of M31 with surface brightness R = 18.2 mag arcsec^. The practical challenges to 
pixel lensing are being overcome by the observers. Theory is not wildly out in front 
of these advances. On the contrary, it is struggling to catch up with them. 



In classical lensing, one does not actually measure A. Rather one measures the 
total fiux F from the lensed "star". 



where A is given by equations (1.1) and (1.3), Fq is the fiux from the star in the 
absence of lensing, and B' is the fiux from any unlensed blended source within 
the star's PSE. If B' is assumed to be zero (as it often is), then Fq can be very 
well determined from the unlensed portion of the light curve. In this case the 
measurement of F does directly yield A. However, since lensing searches are always 
carried out in crowded fields (and since in any event the lens itself can contribute 
unlensed light) the blend term cannot be assumed to be 0. The only quantity that 
is really known from the unlensed light curve is B = B' -\- Fq. Rewriting equation 



2. Pixel Lensing Light Curves 



F{t; to, /3, CO, Fo, B) = FoA{t; to, /3, to) + 5', 
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(2.1) in terms of B gives, 



F{t;to,(3,Lo,Fo,B) = Fo[A(t; to, /3, cu) -1] + B. (2.2) 

Since B is very well known, the observed light curve when written in this form 
must be fit to four unknown parameters, Fo,^o,/3, and io. 

The basic idea of pixel lensing is to subtract from the current image I[t) a 
reference image R derived from images taken before or after the event occurs. The 
actual construction of a reference image is a rather subtle question which will be 
addressed in § 9.6. For the present, I will assume that R is identical to I[t) except 
that 1) the star is not lensed and 2) R has much less noise than I{t). The difference 
image Dit), 

D{t) = I{t) - R, (2.3) 

is then a perfectly blank field except for 1) a PSF with total fiux F = Fo[A — 1) 
and 2) photon noise. That is, the fiux is given by equation (2.2) with B = 0. Since 
i? is a known constant for both pixel lensing and classical lensing, the structures 
of the light curves are essentially the same for the two cases. 

Pixel lensing differs from classical lensing in only one fundamental way. In 
pixel lensing there are many stars per resolution element, so the photon noise is 
dominated by light from stars not being lensed. This means that the noise is 
virtually independent of the magnification. In classical lensing by contrast the 
photon noise is generally dominated by light from the lensed star. Thus, while the 
form of the lensing signal given by equation (2.2) is the same, the form of the noise 
is different. 

For spike pixel lensing there is a second difference. In this regime there are 
so many unresolved stars per resolution element, one is generally able to observe 
events only when /3 ^ 1. The excess fiux then takes the limiting form, 

Fo[A{t;toJ,uj) - 1] ^ ^G'(t;to,cUeff), (2.4) 

P 

where 

and cUgff is the effective time scale. For reasons which will become clear below, I 
refer to G as a "filter function" . Notice that whereas four parameters are required 
to describe an ordinary lensing event, only three (Fq//?, ^0, and coeff) are required 
to describe events in the spike regime. That is, to the extent that the light curve 
is adequately described by G, it is impossible to measure either the unlensed fiux 
Fq or, more importantly, the time scale lo~^ . 
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Figure 1. Three functions which help characterize the relation of the flux excess function \A{t \ to, /?, i^) — 1] 
[see eq. (1.1)] to its limiting form G{t;ti:j,uj / [3) [see eq. (2.5)]. The function ?y(/3) is the correlation between 
(j4 — 1) and G. The fact that ry ~ 1 shows that the G's are excellent filters to search for lensing events. 
However, since G contains no information about the lensing time scale , this fact also means that it 
is difficult to recover to from {A — 1). The function A(/3) quantifies this difficulty: if the event is detected 
with signal to noise Q, then the fractional error in the time scale is Q/K. The function [C(/3)]^''^ is the 
ratio of the true signal-to-noise for {A — 1) compared to what it would be if the excess fiux function were 
actually G. 



Thus, the first step in analyzing pixel lensing is to quantify the difference 
between (A — f ) and G as a function of fi. To this end I define three functions, 
nifi), A(/3), and ({fi). 
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Figure 1 shows the correlation between [A — 1) and G, 



^ / dt{A{t; to, /3, u) - l)G{t; tp, cueff ) 

{/Jt[A(t;to,/3,cu)-l]2}i/2[/^^tG(t;to,cUeff)2]i/2- ^ " > 

Note that > 0.87 over the entire range < /3 < 1. This shows that even 

classical lensing events are reasonably well represented by the filter function (which 
contains no information about the event time scale It implies that even in 

classical lensing, the time scale is a relatively degenerate parameter. It also implies 
that the two-parameter functions G are reasonably efficient filters to use to search 
for microlensing light curves, even when the events have relatively large impact 
parameters. See § 8. 

I define A(/3) so as to quantify the precision of the time-scale measurement. 
Suppose a series of measurements F[ti) are made at uniform intervals which are 
short compared to for a duration long compared to . And suppose that 
each measurement has the same noise 6F. Then the total signal to noise of the 
event Q is given by = (Fo/SF)'^ ~ 1]^- The same measurements can 

also also be used to determine cu, and this determination will have some error Sco. 
I define A = QSlo/uj. That is, A tells us how much more difficult it is to measure a 
time scale than to merely detect an event. Figure I displays A/20. The important 
point is that for 0.05 ^ /3 ^ 0.5, A is approximately constant at ~ 10. This means 
that an event can be detected with good confidence [Q ~ 10) and still yield no 
information about the time scale. Just to make a measurement of the time scale 
at ~ 20% accuracy requires Q ^ 50. 

The last quantity displayed in Figure I, 

^ Jdt[Ait;to,^,uj)-l]' 

/3-2/JtG'(t;to,cUeff)2 ' ^ ' 

is useful for estimating Q. As I will discuss in § 5, for /3 ^ I, oc Fq j fiui. The 
function ^ is a correction factor to this equation for finite fi: oc Fq({(3) / (3lo. 
For fixed Fq and cu, one then finds Slo/uj oc A(/3)[/3/C(/3)]^/l This combinat ion IS 
monotonic in /3 which means that if the time scale can be measured to a given 
precision for a given /3 then the measurement improves at smaller /3. 
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3. Beyond the Crowding Limit 



Let us suppose that a candidate lensing event is identified and that its center 
z = {zi,Z2) on the difference image Dijit) is Icnown. Here denote the integer 

pixei position, whiie z is a pair of reai numbers. Let V'z;zj be the PSF centered 
at z and normaiized so that V'z;zj = 1- Let the mean number of photons per 
pixei in the originai image Ii-j{t) in the neighborhood of z be n. Then the best 
estimate of the number of photons in the stellar image is Opgf V'z;zj-Dij (t) and 
the variance of the estimate is nOpgf , where the resolution element Opgf is defined 

by, 

%li^Y.^hr (3-1) 

For a Gaussian PSF, Opgf = ^ira^ where a is the Gaussian half-width. For Nyquist 
or "critically sampled" data (cr = 1 pixel), Opgf ~ 12.6 pixels. (I use the same 
symbol Opgf to express the resolution element in pixels or arcsec^.) 

Note that Opgf can vary over the image on two distinct scales. First, because 
of discrete pixelization, the shape of the PSF will depend on where the center is 
relative to the pixel structure. If the data are oversampled, then the pixelized PSF 
will vary smoothly and the sum in equation (3.1) will be well approximated by an 
integral. In this case, Opgf varies very little with position. The more undersampled 
the data, the more Opgf will vary. However, even for the significantly undersampled 
Hubble Space Telescope (HST) WFC2, for which I find a mean (Opgf^ = 8.5 pixels, 
the standard deviation around this mean is only 11%. For purposes of this paper, 
such small variations can be ignored. The resolution element can also vary across 
the field because of telescope optics. However, since all analyses presented here 
will be conducted locally, this variation can also be ignored. 

If the density of detectable stars on an image is 0"^^, it is said to be "crowding 
limited". Below the crowding limit, an increase of the exposure time increases the 
number of resolvable stars as more of the luminosity function rises above the noise. 
Once this limit is reached, however, further increases in exposure time has only 
a marginal impact on the number of resolvable stars. In classical lensing studies, 
the maximum number of stars which can be monitored is equal to the number of 
resolution elements, A^res, 

iVres = (3.2) 
' 'psf 

where Occd is the area of the CCD. Classical lensing is therefore optimized if the 
sampling is just critical (to maximize A^res) and if the exposure times are just 
sufficient to reach the crowding limit. 
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The optimal criteria for classical lensing are not determined by the nature of 
the lensing events, but by the nature of the data reduction programs that are 
employed to detect them. These programs have their origin in the problem of 
doing photometry in crowded fields, such as globular clusters. The goal of such 
observations is to measure the fiuxes from as many stars as possible. However, 
measuring the fiux from individual stars is completely irrelevant in microlensing 
studies. All of the information is contained in the fiux differences Fo[A — 1). 
These italicized claims seem to contradict common sense because there certainly 
would be additional information if one could measure the unlensed fiux Fq. As 
I have emphasized, however, what is measurable in crowded fields is not Fq but 
B = Fo + B'. See equations (2.1) and (2.2). 

The "crowding limit" is not a fundamental barrier to pixel lensing. Deeper 
exposures increase the number of stars per resolution element that are effectively 
monitored. In principle, this number is limited only by the base of the LF at 
Mj ~ 14. In practice, the value of obtaining deeper exposures saturates at a 
certain point. This point is a function of the characteristics of the galaxy (distance, 
surface brightness, LF) and also of the telescope. For some galaxies (like the bulge 
of the Milky Way) this point of diminishing returns is in fact the crowding limit. 
For more distant galaxies, it is well beyond the crowding limit. 



4. The Fluctuation Magnitude Star 

There are two fundamental parameters which characterize pixel lensing. The 
first is the resolution element Opgf given by equation (3.1). The second is the 
"fiuctuation fiux" F^, the ratio of the first two moments of the LF <f>{F)^ 

_ JdFHF)F' 
*- JdF<P{F)F- ^ ' 

The fiuctuation fiux (and the corresponding fiuctuation magnitude) are of course 
also fundamental to the surface-brightness-fiuctuations (SBF) method of measur- 
ing distances (Tonry 1991). For / band, the absolute fiuctuation magnitude is 
Mj = —4.84 + 3.0(y — /). The distance modulus is determined by measuring the 
fiuctuation fiux, converting to an apparent fiuctuation magnitude, m/, correcting 
for extinction A/, and then subtracting fi = rhj — Mj — Aj. In lensing only the 
empirically determined m/ comes into play. 

The fiuctuation magnitude star enters almost every aspect of pixel lensing. As 
I show in the next section, the detectable event rate in the spike regime is directly 
proportional to the rate of photon detection from a fiuctuation magnitude star. 
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In the following section I show that helps define the boundary between semi- 
classical and spike regimes. In § 9, I show that the relationships between photon 
noise and various types of systematic error can be quantified in terms of F^. 



5. Spike Pixel Lensing 

Consider a galaxy with surface brightness S observed when the sky brightness 
is 5'sky. Define /€ = (1 + Ss]^y / S)~^ . Let a be the number of photons collected by 
the telescope per unit time per unit fiux, and let t^x be the exposure time. The 
signal to noise of a measurement of excess fiux due to a lensing event at time is 
then, 

g. = Fo[A(t.)-l](^)'^'. (5.1) 

(I have ignored readout noise here. However, it could easily be included in this 
formalism by augmenting the sky noise. For example, if the readout noise is 10 e~ 
per pixel, the sky counts could be increased by 10^.) 

The total signal to noise of the event is then given hy Q = Qf)^^^ ■ I 
now make three simplifying assumptions. First, I assume that the observations 
are carried out uniformly at intervals small compared to ■ Second, I assume 
that the sky noise and hence k, is constant over time. Third, I assume that Opgf is 
constant. The sum can then be converted into an integral which yields. 



irF^aK C(/3) 



(5.2) 



where a is the photon detection rate averaged over the duty cycle of the telescope 
and ( is the suppression function given by equation (2.7) and shown in Figure 1. 

The total signal to noise is a monotonic function of /3, so that for fixed Fq and 
cu, and for a given threshold of detectability Qmin, there is a maximum impact 
parameter /3max such that events below this value are detectable and events above 
it are not. Note from Figure 1 that ^ — )• 1 in the limit /3 ^ 1 and that ( falls 
rapidly for /3 ^ 1/4. For simplicity I approximate ({/3) ~ 0(1/4 — /3), where is 
a step function. The maximum impact parameter is then, 

/3^,,(Fo) = mini ,U. (5.3) 



' mm 



^cuS'Opsf 4 



The rate of detectable events from a given star is just /JmaxTo where Fq is the 
classical event rate given by equation (1.5). Integrating this rate over the LF and 
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the area of the CCD, I find 



r = ■^TN.esnT^, r, = aF^ (5.4) 

^min 

where is the photon detection rate from a fiuctuation magnitude star, and ^ is 
a correction factor arising from the the upper hmit /3niax < 1/4, 

Here -Fmin is the minimum fiux hmit which results from finite-source size effects. I 
will discuss these effects below in § 7 but for the moment ignore them. 

The fiux -Fmax represents the boundary between two regimes. Lensing events 
of stars with Fq > -Fmax can be recognized even when they have low impact pa- 
rameters (/3 ^ 1/4)- The light curves associated with these events probe the entire 
Einstein ring and therefore tend to give good information about the timescale lo~^ 
and the stellar fiux Fq. Stars with fiux Fq < -Fmax give rise to events that are 
detectable only as spikes near the center of the Einstein ring, and which therefore 
do not generally give information about the timescale. 

If -Pmax is brighter than -Fum, the brightest star in the LE, then ~ f and pixel 
lensing is in the spike regime. The event rate is then related to r through directly 
observable quantities, F^, A^res, and Qmin- This means that the optical depth 
can be determined from the observations. Note that in classical lensing the optical 
depth determination requires measurements of uj~^ [see eq. (1.4)]. By contrast, in 
the spike regime of pixel lensing the optical depth is determined with essentially no 
information about the individual time scales of the events. Since /3 < f /4, the light 
curve is well described by the filter functions (^(t; to, cUgff ) which does not depend 
on id. See r][l3) in Figure f . 
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6. Semi-Classical Pixel Lensing 



The boundary between the spike and semi-classical regimes is the point where 
= i^max- Exphcitly, 

rV-^gk, N.^'^. (6.1) 

Note that the left hand side is the ratio of the number of photons collected from 
a fluctuation magnitude star during an Einstein crossing time to the number of 
fluctuation magnitude stars in a resolution element. If the exposures are much 
less deep than this, then ~ 1, the case discussed in the previous section. If 
the exposures are much deeper than this (the semi-classical regime) then ^ 1. 
Equation (5.4) remains formally valid, but is no longer a very useful way to express 
the event rate. 

In particular, what one gains from deeper exposures is different in the semi- 
classical and spike regimes. Eor spike pixel lensing, the total number of events 
rises directly with a/Opgf. Eor semi-classical lensing, the event rate grows much 
more slowly. This is well illustrated by the simple quasi-realistic LE (f>{F) oc 
F-'Q{F]i^ - F), for which ( = {F^.jF]i^f[l + 2 ln(Fii^/Fniax)], implying F oc 
[1 -|- 2 ln(Fiiin/Fniax)]- The growth of F with a/Opgf therefore goes from linear for 
Fmax > -Plim to Square root at -Fmax ~ -Plim to logarithmic at -Fmax ~ F^^ (typically 
~ 1 mag below -Fiim)- Hence, from the standpoint of detecting as many events as 
possible, it would be better to use the available telescope time to look at another 
galaxy or perhaps another field in the same galaxy. Recall from § 2, however, that 
if the exposures are sufficiently deep, one can measure the time scale of individual 
events. The time scales can be used not only to determine r, but to measure the 
mass spectrum of the lenses as well (Han & Gould 1996). One can also measure 
the individual fiuxes Fq of the source stars. Microlensing is probably not the most 
efficient method for studying the EE of the nearest galaxies, but may be the only 
method for distant ones. Hence the source fiux measurements could prove to be 
uniquely valuable. 

The minimum signal to noise required to measure the fiux with ~ 20% accuracy 
is Q ^ 50. As I will show in § 8, the minimum threshold for mere detection is 
Qmin ~ 7. Thus, in order to measure time scales for even the brightest source- 
star events, one must reach two mag below the tip of the LE, -Fmax ~ 0.15Fiini. 
Typically this is ~ 1 mag below F^. Hence, there is a roughly 2 mag interval 
between the purely spike and purely semi-classical regimes. Equation (6.1) defines 
the midpoint of this transition. 
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While a 2 mag interval seems modest by astropliysical standards, moving 
through this interval requires an improvement of a factor ~ 50 in telescope time, 
telescope size, or PSF size. To a first approximation, then, whether pixel lensing 
is in the semi-classical or spike regime depends primarily on the galaxy's distance 
and only secondarily on the observing program. 

7. Finite Size Effects 

The finite angular radius of the star 6s fundamentally limits pixel lensing. If 
/3max < ^s/2^e, then the maximum magnification will be underestimated by the 
point source formula (1.1), and the event will not in fact be seen (Gould 1995b). 
This occurs when F < -Fmin where, 

^min = '^'^tsF^^^, (7.1) 

and ts = 6s/io6e is the crossing time of the star. Note that equation (7.1) is actually 
a self-consistent condition because the left hand side contains the fiux of the star 
and the right hand side contains its radius. 

One can draw a few general conclusions from equation (7.1). First, if one is 
attempting to detect lenses in or near the source galaxy, then ts ~ Vs/v where 
Vs is the source radius. Near the tip of the LF, ~ 0.5 AU. For typical stellar 
mass ~ O.3Af0 and typical galaxy or cluster internal distance d\s ~ 1-100 kpc, the 
physical Finstein radius is c^ol^e ~ 2-20 AU. That is, i^maxZ-Pmin ~ 2-6. Hence, 
for spike pixel lensing to be observable at all, -Fmax must be no more than one or 
two mag brighter than -Fiim- Recall from § 6, that if Fjim is more than one or two 
mag fainter than -Fum, then the pixel lensing is semi-classical. This implies that 
total range of spike pixel lensing is only a few mag. An important exception to 
this rule occurs for Galactic lenses. In this case Be is typically larger by a factor 
dos{d\sRo)~^^^ where i?o is the galactocentric distance, so that -Fmin is almost always 
negligible. This implies that it is possible to use external galaxies as sources to 
search for Galactic lenses, even when lenses within those galaxies are undetectable. 

If the time interval between exposures is longer than t^, then this interval enters 
equation (7.1) in place of ts- In fact, for ground-based observations from a single 
site, the situation is considerably more complicated. Observations during a single 
night would be sensitive to events with effective time scales uj~^ ^ 4hrs, while 
observations over several days would be sensitive to those with uj~^ ^ 1 day, but 
the sensitivity would be poor for intervening time scales. 
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8. Searching For Events 



Suppose that the images all have the same PSF and all have the same 5'sky. 
One might then imagine searching for events by forming sums of the difference 
images Diti), 

D{h,Lo,^) = Y,D{U)G{U-h,Lo,^). (8.f) 

i 

For an event with /3 ^ f , effective time scale exactly equal to cUgff , and peak exactly 
at to, such an image would contain a PSF with signal to noise equal to the Q of the 
event itself. One could search the image for point sources and all such events would 
appear. Fven for finite /3, the signal to noise would only be reduced by 
[see Fig. f and eq. (2.6)], a modest correction at most. If the peak differed from to 
by a time ^ lo~^ or the time scale differed from lo~^ by a small fractional amount, 
the signal to noise would also fall by only a small amount. 

This search method is a practical one for some observing programs, especially 
those conducted from space. In any event, the method can be used to estimate 
the number of independent event searches which are being made. The overlap 
integral between two filters (^(t; to, cUgff ) and (^(t; t'g, cu^^) will be high provided 
l^eff (^0 ~ ^o)| ^ 1 and | ln(c<;eff /cu^^) | ^ I. Hence there will be ~ A^obs In^obs inde- 
pendent combined images _D(to, cUgff ), where A^obs is the ratio of the duration of the 
observations to the shortest effective time scale being searched. Note that A^obs will 
typically be equal to the number of observations. The number of independent reso- 
lution elements on each image is A^res- Hence, a total of ~ A^resA^obs In ^obs indepen- 
dent searches are conducted. For an ambitious program, one might have A^res ~ 10^ 
(for moderately oversampled data on a 2048^ chip) and A^obs ~ lyr/Ihr ~ 10^. 
This implies ~ 10"^"^ independent searches. The threshold required to avoid detec- 
tion of noise-induced events is then Qmin ~ 7. 

A more generally useful search technique would be to form images D' iti) by 
replacing each pixel with the value obtained by fitting a PSF centered at that pixel 
to the image. Then for each pixel (j, k) and each filter function G{t] to,c<;eff), form 
the signal-to-noise ratio. 



^ 



psf ,i 



-^1'' D\^{ti)G{U■U.^^^y^ 



^ 



f8.2l 



psf ,i 



where n^]^ is the number of photons per image in pixel (j, due to the galaxy 
(ignoring sky). Finally search for pixels (j, A;) and filters at (to,c<;eff) for which 
Q ^ Qmin- In practice, one would actually set the threshold slightly below Qmin 
and then search for the local peak in Q in the neighborhood of the discrete peak 
at (j, A;, to, cUgff ). It is this peak which should exceed Qmin- 
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9. Noise From Systematic Effects 



So far I have simply assumed that the only obstacle to detecting and measuring 
pixel-lensing events is photon noise, an assumption that some might politely label 
"optimistic". For example, I assumed that the images can be aligned perfectly 
with no loss of information and, perhaps more strikingly, that the PSF in the 
reference image R is identical to the one in the current image I{t). I assumed that 
R contains no noise and that cosmic ray events (CRs) play no role. Simply to list 
these complicating factors is enough to convince most people that pixel lensing is 
hopeless. However, it is possible to evaluate each effect quantitatively and to show 
that they do not compete with photon noise for realistic observing programs. 

9.1. Basic Formalism 

Consider a dense star field composed of stars with angular positions Zj and 
fiuxes F^. The surface brightness as a function of position is then given by. 



As seen from Farth, the surface brightness will not appear as the discrete sources 
represented by the (5-functions in equation (9.1), but as smeared out by the PSF 
t/j with a full width half maximum (FWHM) ^see which is generally of order an 
arcsec. The apparent surface brightness is then given by the convolution. 



If the typical separation of the stars is much smaller than ^see, then the stars will 
be unresolved. This is the regime of pixel lensing. 

Let us suppose that the field is imaged twice, first at time ti when none of 
the stars are microlensed (the reference image i?), and later at time ^2, when 
the jth star is microlensed with magnification A. Assume for the moment that 
one can align the two images photometrically and geometrically (see below). The 
first image is subtracted from the second, and the difference image has surface 
brightness AS*, 



A5'(z) = S{z, t2) - S'(z, h) = - i^i)oSo + {A- l)F,ij2oS{z - z,), (9.3) 



where t/ji and t/j2 are the PSFs at the two epochs. If we further suppose that 
V'l = V'2, then the first term in equation (9.3) disappears. The difference image 
is just an isolated point source convolved with the PSF. This is the situation I 





S = tpoSo. 
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assumed in § 3 where I evaluated the photon noise -F-y noise (assuming that R is 
noiseless), 

Recall that t^x is the exposure time, a is the photon detection rate, and k, is the 
ratio of flux generated by the galaxy itself to total flux. Note that I have redefined 
Opsf as a continuous integral rather than a discrete sum. I consider the effects of 
discrete pixelization in § 9.5 below. 

9.2. Noise From Time Variation of the PSF 

First, I compare the noise due to time variation of the PSF with the photon 
noise -F-ynoise- If the two PSFs in equation (9.3) are not the same, then the first 
term AS'^''^ = {1^2 — V'l)°'5'o will be "bumpy". 

Suppose that the difference image contains a stellar profile at a given position, 
(0,0), and suppose that the fiux in the star is measured by PSF fitting, i.e., by 
multiplying the PSF by the observed fiux. Then the residual bumps AS'^''^ will 
add noise to this measurement. To evaluate this noise, first consider a patch of the 
image at position z and with angular area AO which is large enough to contain 
many stars, but is still small compared to the PSF. If the patch contains an excess 
fiux F, then the PSF fitting procedure will mis-estimate the fiux of the star at 
(0, 0), by an amount F {At/j^li^) / (V'lV')? where 

AV'z(z') = V'2(z'-z)-V'i(z'-z), (9.5) 
and where the Dirac brackets indicate an inner product, 

(i^lV') = J SzH{z)i^{z). (9.6) 

Note from this definition and equation (9.4) that Opgf = (V'lV') ^ ■ L^t us assume for 
the moment that the mean surface brightness S of the galaxy is generated entirely 
by fiuctuation magnitude stars with fiux F^. There will be a total of SAil/F^ such 
stars in the patch, and hence the square of the error induced by the patch will be 

A(i^p'sfnoise) = ^[F^^vsi (A^'zlV')]'- (9-7) 

In fact, one may show that this formula remains valid for an arbitrary luminosity 
function. Integrating equation (9.7) over the entire PSF and dividing by equation 
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(9.4) yields 




psf noise 



'psf 



/ 



^(AV'zlV')', 



(9.8) 



where is the number of photons collected during the exposure from a fluctuation 
magnitude star. For two Gaussian PSFs with standard deviations ai and (72, and 
(5(7 = (72 — (71 < (7, Opsf / (fz (AV'zlV')^ ~ ('^^/^)^/4, SO that 



In a f hour /-band exposure of M3I, aim telescope would collect ~ 1500 
photons from a fluctuation magnitude star (m/ ~ 23.5). Hence the PSFs must 
be similar to within ~ 5% if the PSF noise is not to dominate. It is possible to 
measure the PSF from foreground stars to much better than this accuracy. Using 
this PSF measurement, the reference image can be convolved to the seeing of the 
current image (see § 9.6). 

What underlying physics permits one to eliminate noise from the PSF differ- 
ence, A?/;? The noise is generated by the SBF "bumps" in the surface brightness 
of the galaxy. To the extent that these bumps are convolved with different PSFs, 
they leave behind a noisy difference image. There are fluctuation magnitude 
stars per resolution element [cf. eq. (6.1)], hence the bumps have flux ~ NJ F^, 
which falls inversely with distance. For galaxies at least as far as MSI, there are 
always field stars that are much brighter than the bumps. The PSF can therefore 
always be measured more accurately than is required to remove the noise from 
bump residuals. 

The difference image consists of an isolated point-source PSF on a fiat back- 
ground (much like the image of a random high-latitude field) and hence it is tempt- 
ing to estimate the errors in the fiux using the same procedure one would use for 
isolated-star photometry: first measure the variance of the background pixel counts 
and then multiply this by the number of pixels in Opgf . The product should then 
be variance of the fiux estimate. However, if one applies this standard approach 
to the problem of PSF-generated noise, one overestimates the variance by a factor 
~ 8. Specifically, one finds analytically that the variance of the difference im- 
age is given by var^S"^'^) = F^S J d"^ z[A'ijj{z)]'^ . If this variance is then measured 
and naively multiplied by Opgf , one obtains a result which is 8 times greater than 
equation (9.9). It is easy to understand the source of the overestimate. In a dif- 
ference between two images with similar but not identical PSFs, stars will have a 




7 noise 



(9.9) 
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characteristic profile, e.g. "rings" or "shadowed mountains". However, because the 
net flux in these structures is zero, and because they are anti-correlated on scales 
of order the PSF, they actually induce smaller errors than would fluctuations of 
the same amplitude but randomly distributed. Noise from geometric alignment, 
photometric alignment and pixelization (see below) are also overestimated by the 
standard procedure. 

In any event, this calculation shows that it is not straight forward to estimate 
the flux errors from the gross noise characteristics of the difference image. Instead, 
error estimates should be made from repeat photometry at short time intervals and 
from the scatter of periodic variable star fluxes around their mean light curves. 

9.3. Geometric Alignment 

The foregoing calculation immediately permits one to estimate the accuracy 
required for the geometric alignment. In this case, t/j2 is the same as t/ji but 
offset by some amount 69. One then finds Opgf J (A?/'z|?/')^ ~ ((5^/(t)^/4, so 
( -^geom noise/ -^7 noise )^ ~ Kn^{S9 / a)"^ / 4. In the above example, this implies a re- 
quirement to reduce 66 ^ 0.05(7 or about I/IO of a pixel for mildly oversampled 
data. In practice, one can do much better. The underlying physical reason is 
the same as above: foreground stars are much brighter than the SBF bumps, so 
residuals in the bumps arising from uncertainty in the image position are small. 

9.4. Photometric Alignment 

After two images Ii and I2 are aligned geometrically, they will still suffer 
some photometric misalignment due to differences in the sky level and atmospheric 
extinction. That is /2 = a + bli. An error 6b in estimating the linear term in this 
transformation can produce noise in the difference image. (Photometry is not 
affected by a constant offset.) The amplitude of this noise will simply be 6b/b 
times the amplitude of the SBF (or "graininess") of the original image. The SBF 
can be evaluated by substituting At/j — )• [6b/b)t/j in the previous analysis. This 
then implies a ratio of alignment noise to photon noise (-F'align noise/ -^7 noise )^ ~ 
Kn^{6b/b)'^ /2, yielding for the above example the easily reached requirement 6b/b ^ 
3%. For galaxies at least as far as M3I, photometric alignment is most easily 
accomplished by rescaling the images so that their means and standard deviations 
are equal. The "signal" for this transformation is the fiux from the entire galaxy 
which of course is much larger than the fiux from any of the bumps that one is 
trying to align (Melchior 1995). 
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9.5. PiXELiZATiON Noise 



Finally, there is noise induced by pixelization. If the images are aligned by 
shifting them using a flux- conserving algorithm (like linear interpolation) and if 
the excess flux [A — l)Fj were measured by aperture photometry, then pixeliza- 
tion would not induce any noise. However, aperture photometry itself introduces 
an enormous amount of noise. PSF fitting minimizes photometric noise and was 
assumed when I analyzed the sensitivity of pixel lensing in § 5. To understand pix- 
elization noise, suppose a star with unit fiux is centered in an image at (z, j) = (0, 0) 
and suppose that another image is shifted relative to the first by (a, /3) fractional 
pixels. The fiux at a point (x, y) in the neighborhood of the center of pixel (z, j ) can 
be expressed as a Taylor series V'(^7 v) = i^ihj) + {x — i)il)x{i,i) + (j/ — j)4'y{h j)--- 
where subscripts represent derivatives. Using this expression, I find that if one uses 
linear interpolation on the second image to approximate the fiux that would have 
fallen in (z,j) in the first, one makes an error, 

AV'(^, j) = -^Kl - «)V'xx(^, j) + /3(1 - f3)i^yyi^, j)], (9.10) 
Within the framework of the continuous-variable formalism already developed, 

AV'z(z') = -^[«(1 - «)V'xx(z' - z) + /3(1 - /3)V',,(z' - z)], (9.11) 

I will use equation (9.11) rather than equation (9.10) primarily because it is more 
convenient. However, I should note that equation (9.11) is actually more accurate 
because it satisfies exactly the fiux conservation condition, J d"^ z' Atpz^z') = 0, 
while equation (9.10) salsifies the analogous discrete equation only approximately. 
Assuming a Gaussian PSF and combining equations (9.8) and (9.11) yields. 



pixel noise KH,^ 



F2 . 128a' 

7 noise 



■{3[a(l - a)f + 2a{l - a)/3(l - /3) + 3[/3(l - (9.12) 



For an ensemble of images, a and /3 will be uniformly distributed over the interval 
(0, 1/2) yielding a mean square error. 



pixel noise KH,^ 



F2 . 500(7^ 

7 noise 



(9.13) 



Fquation (9.13) implies that there is a strong incentive for obtaining highly over- 
sampled data. An alternative would be to reduce by obtaining many shorter ex- 
posures. This will be necessary in any event in many cases to avoid over-exposure. 
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I should note that in cases ol excellent pointing (e.g. HST wheie a, /3 ^ 0.05 pixels) 
equation (9.13) should be multiplied by ~ 30 (ct^)- Equation (9.13) has important 
implications for follow up observations of classical lensing events which I discuss 
in § 10.5, below. 

From the discussion at the end of § 9.2, the ratio of true noise generated by 
systematic effects to the naive estimate one would make based on the variance in 
the pixel counts is given by 



For the PSF variation, geometric misalignment, photometric misalignment, and 
pixelization (and for Gaussian PSFs) this ratio has values of 1/8, 1/2, 1/2, and 
1/8, respectively. Thus, all systematic effects have less impact on the photometry 
than would appear based on the measured variance in pixel counts. 

9.6. Reference Image Construction 

I have so far assumed that the noise in the reference image is negligible com- 
pared to the noise in the current image. From the standpoint of photon statistics 
alone, this implies that there must be many more exposures which are used to 
construct the reference image than are used to track the event. For events that are 
of order the duration of the observations this is obviously impossible. The noise 
level for these long events will be augmented by a factor ~ 2 relative to what I 
have assumed. However, for all but the shortest observation programs, the great 
majority of events will be short compared to the duration of the experiment, so 
that photon noise in R will not generally be a limiting factor. 

However, it is important to recognize that the requirements on R are not triv- 
ially satisfied. For example, it would not be adequate to form a single reference im- 
age out e.g., 5 "good" images that are not affected by the event, and then subtract 
this from each event image. In this case, the difference images Diti) = I{ti) — R 
would individually be dominated by the photon noise of /(tj), but in the convolved 
image used to search for events _D(to, cUgff ) = G{ti] tQ,LOf.f[)D{ti), the noise from 
R would still dominate. This is because the contributions from G{ti)I{ti) would 
add incoherently, while the contribution from G{ti)R = R^^ G{ti) would add 
coherently. Fither each image I{ti) must have its own reference image constructed 
of several images not in the event and not used to reference other images, or a 
single reference image must be constructed from a very large number of images. 

These requirements lead to two different types of solutions. The first is to group 
the reference images by seeing and pair each group with a single I{ti) which has 




(9.14) 
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somewhat worse seeing than the group. Then convolve each of the members of the 
group to the seeing of I{ti) and form their mean excluding 3 a outliers. (One could 
also form the median, but this has (it/ 2)^^"^ more noise without, in my view, any 
compensating advantages.) A second solution, applicable mainly to space-based 
observations where the seeing is constant, is to form a single reference image from 
all the images. In either case, the noise of the reference images can be reduced well 
below that of the /(tj). 

CRs are automatically removed from R by the above procedure. A more del- 
icate problem is to remove them from the /(tj). The problem only arises when 
searching for events for which cOefftex <^ 1- For longer events one can remove CRs 
in the standard way by comparing images at neighboring times. Even for the 
shorter events, CRs pose no problem once a candidate event has been identified: it 
is easy enough to check whether a particular measurement has been corrupted by a 
CR. CRs are only a problem in the search for short-event candidates. By searching 
the Diti) for spatially correlated 3(7 outliers whose global structure is inconsistent 
with a PSF, one could probably remove essentially all CRs. However, it may be 
that CRs become a limiting factor for short events in some cases, particularly for 
long, space-based exposures. 



10. Applications 

The analysis presented here can be used to evaluate prospective pixel lensing 
experiments and to develop new pixel lensing ideas. Here I give a few examples. 

10.1. Classification of Target Galaxies in / Band 

For simplicity I assume the fiuctuation magnitude is Mj ~ —1.2 in all galaxies, 
and that a telescope with diameter D collects 10(-D/1 m)^ photons per second from 
an / = 20 star. Critical pixel lensing [F^ = -Fmax) can be evaluated from [cf. eq. 
(6.1)], 



i^max VMpcy ImVhry V 1" 7 V 7 

where to^s is the observation time per day, and S'lg is the total surface brightness 
(including sky) in units of 19 mag arcsec""^. 

Thus dedicated observations with a 2 m telescope would bring the disk and all 
but the innermost part of the bulge of M31 well within the semi-classical regime. 
If a substantial fraction of the halo of M31 is composed of massive compact halo 



21 



objects (MACHOs), it should be possible to measure not only their total optical 
depth but the individual time scales of many of the events. 

Even if there are no MACHOs, there will be lensing of giants in the M31 bulge 
by the more numerous dwarf stars in the bulge. Because the bulge of M31 is semi- 
classical, it is possible to measure not just the time scale of these events, but also 
the unlensed flux from the giant sources. That is, M3f bulge-bulge lensing can 
yield information on both the mass function and the LF of the M3f bulge. Han 
(f996) has applied the present formalism to analyze pixel lensing of M3f . 

There are several other large galaxies within a few Mpc, such as M81 and 
Cen A. Dedicated observing programs on 4 m class telescopes with excellent seeing 
(^see ~ 0."5) could bring these into the semi-classical regime and would allow study 
of their mass and luminosity functions. They are accessible in the spike regime with 
more modest programs which would allow measurement of their optical depths. 

As I have elsewhere discussed (Gould 1995b), it would be particularly interest- 
ing to search for pixel lensing toward M87 because of the possibility of detecting 
intra-cluster MACHOs in the Virgo cluster. Dedicated observations with iJ^T just 
barely bring the outer part of M87 into the semi-classical regime. Most of the 
galaxy would be in the spike regime and the observations would yield information 
about the total optical depth, but not about individual time scales nor about the 
LF of M87. 

Coma is much richer than Virgo and would therefore be an even more inter- 
esting cluster to probe for intra-cluster MACHOs. However, the scaling of 
equation (10.1) together with the fact that Coma is ~ 5 times farther than Virgo 
makes such a search impossible with existing equipment. 

The Milky Way bulge, and the Magellanic Clouds are all so close that even 
with the 2-10 minute exposures of the current searches, they are well within the 
semi-classical regime. This is of course to be expected since the current programs 
are searching for classical lensing. 

10.2. Space-Based versus Ground-Based Pixel Lensing 

While it is not possible with existing equipment to survey galaxies much beyond 
M87 for pixel lensing, such studies may be possible in the future. For diffraction 
limited images obtained from space, the distance Jqs at which pixel lensing can 
be observed scales directly with telescope diameter, D. To see this, note that 
L^ oc (D/dos)^ and Opgf oc . For fixed number of pixels in the detector (or 
fixed physical size of the galaxy in the case that the field is larger than the galaxy) 
Occd oc . Equations (5.4), (6.1), and (7.1) are all invariant under these changes. 
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On the other hand, for ground-based telescopes of fixed PSF, the improvement 
with larger aperture is much more gradual. For galaxies contained within a single 
CCD frame, A^res oc . Hence the event rate [eq. (5.4)] in the spike regime and 
the boundary between the two regimes [eq. (6.1)] stay constant if the aperture 
changes with distance as D oc (compared to D oc Jqs from space). While the 
comparison is somewhat difficult, I estimate that for pixel lensing of M87 dedicated 
observations on Keck would be ~ 13% as effective as dedicated observations on 
HST. Pixel lensing of distant galaxies therefore appears to require space-based 
observations. 

Various techniques are now being developed to improve resolution from the 
ground including adaptive optics and optical interferometry. These techniques are 
probably not useful for pixel lensing which requires a wide field of view over which 
the PSF can be measured locally with good precision. 

10.3. Infrared versus Optical Pixel Lensing 

Based on a Baade's Window K band LF provided by G. Tiede (private com- 
munication, 1995), I estimate a fiuctuation magnitude ~ —5.7, i.e., ~ 4.5 mag 
brighter than in /. I estimate Afjy ~ —5.4. Ignoring relative detector efficiencies, 
this implies that F^ is ~ 16 times higher in K than /, and 20 times higher in 
H . However, this substantial advantage of infrared (IR) over optical pixel lensing 
is undercut by a number of disadvantages. First, for ground-base observations, 
galaxies are well below the sky in K and i7, but are mostly above the sky in /. 
Second, in the IR it is usually necessary to chop between the object and the sky 
in order to properly subtract the variable sky. This may not be essential for rela- 
tive photometry in K band, but time- variable excitation of water lines in H leads 
to fringes which shift their angular position with time unless sky frames taken at 
similar times are subtracted. While the angular scales of these fringes are many 
arcsec and hence much larger than the PSF, the amplitudes are of order the galaxy 
surface brightness. Unremoved fringes would therefore make pixel lensing impossi- 
ble. Since IR arrays suffer from a number of other problems including smaller size 
and greater non-linear response than CCDs, ground-based IR pixel lensing does 
not appear promising. 

However, for space-based observations, H band would have substantial advan- 
tages over / band assuming similar quality detectors were available. From equation 
(5.4), F oc F^/c/Opsf. For diffraction limited PSFs, (F^/Opgf )jy/(F^/Opsf )/ ~ 5. In 
the inner parts of galaxies /€ ~ 1 for both / and H . In the outer parts, the sky in / 
cannot be neglected, particularly for Virgo and Coma which are near the ecliptic 
and hence prone to scattered zodiacal light. Since this scattering is weaker in i7, 
the IR also has a higher n. IR observations from space may therefore be the best 
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approach to pixel lensing in the future assuming that present-day problems with 
detectors are overcome. 

10.4. Variable Stars 

Variable stars are the most important background for classical lensing searches, 
and this is also true of pixel lensing. For spike pixel lensing, variables pose only a 
minor problem because the events have flux changes corresponding to Mj ^ —4.5 
with characteristic times uj~^ ^ 4 days. Variables with these characteristics are 
rare. However, as pixel lensing moves into the semi-classical regime, it becomes 
sensitive to lower-flux events and hence to contamination by a wider class of vari- 
ables. On the other hand, in this regime the light curves are more accurately 
determined so there is a greater basis for discriminating lensing events with their 
characteristic profile given by equation (2.2) from variables stars. Nevertheless, 
the only way to take systematic account of the variable background is to measure 
the frequency of various classes of variables. This can be accomplished partly by 
searching for variables in the galaxy, a search which grows naturally out the lensing 
search itself. Additional information comes from variable searches in more nearby 
galaxies (which can be done as part of other lensing studies). The higher signal-to- 
noise ratio for the same luminosity star then allows one to unambiguously identify 
a variable in a nearby galaxy which would be difficult or impossible to pick out in 
a more distant galaxy. Variable searches in the Galactic bulge thus serve as a basis 
for estimating the types and frequencies of variable contaminants of M31 lensing 
studies. M31 could play the same role for M87. 

In addition, variables are of interest in their own right. All four classical lensing 
experiments have produced extensive variable catalogs, only a small fraction of 
which have been published (Udalski et al. 1994b, 1995; Alcock et al. 1995d; Grison 
et al. 1995; AUard 1995). 

Pixel lensing studies can also be used as sensitive searches for variables. As 
an example consider a 3-day Cepheid with mean mag / = 21.5 in the inner disk 
of M31 at a projected position where the surface brightness of the M31 bulge is 
/ = 19 mag arcsec"^. Suppose that a pixel lensing study were carried out with 1 
hour observation per night on a 1 m telescope in l" seeing. The mean mag of the 
Cepheid is ~ 6 times brighter than Af/, but since there are ~ 200 fiuctuation 
magnitude stars per resolution element, the Cepheid has only ~ 1/2 as much 
fiux as a typical SBF. If the Cepheid were re-observed when its fiux had changed 
by ~ 30%, the difference would amount to ~ 15% of a typical SBF, and would 
be extremely difficult to recognize by standard techniques. Indeed, no such short- 
period Cepheids have been discovered within the bulge of M31. However, on a pixel 
lensing difference image, this fiux change would yield a ~ 5 cr PSF, which could be 
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detected. More important, by searching the whole series of images with a variable- 
star filter, the Cepheid would show up very strongly. It would not be necessary 
to search with a filter exactly matched to the Cepheid light curve. Sinusoidal 
filters could easily pick out most periodic variables. After they were identified 
as variables, the forms of their light curves could be more precisely measured. 
The same method could be used to detect much fainter variables. In fact with a 
somewhat more aggressive program (4 hours per night on a 2 m telescope) one 
could hope to detect RR Lyraes in the bulge of M31 over the course of a season. 

10.5. Pixel Methods Applied to Classical Lensing Follow Up 

At present two groups have organized follow-up observations of events found in 
classical lensing searches (Pratt 1995; Sackett 1995). By obtaining more accurate 
photometry with shorter sampling times, they hope to detect planetary or binary 
companions to the principal lens (Mao & Paczyhski 1991; Gould & Loeb 1992; 
Udalski et al. 1995b; Bennett et al. 1995) measure the proper motion of the lens 
(Gould 1994; Nemiroff k Wickramasinghe 1994; Witt k Mao 1994), and possibly 
its parallax as well (Gould 1992; Alcock et al. 1995c). These measurements often 
require, or at least would benefit from, photometry with ^1% accuracy which is 
generally considered the limit for standard techniques in crowded fields. 

In contrast to the classical searches themselves, the follow-up observations do 
not need to cover large areas and therefore suffer no pressure toward Nyquist 
sampling. If pixel-lensing techniques (including linear interpolation are employed) 
then the limit on accuracy due to pixelization is given by equation (9.13), with 
replaced by the number of photons collected from the lensed star. Thus, these 
searches should be conducted with at least 5 pixels per FWHM (cr ^ 2 pixels). 
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